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The role of the tip in inelastic electron tunneling spectroscopy (lETS) performed with scanning 
tunneling microscopes (STM) is theoretically addressed via first-principles simulations of vibrational 
spectra of single carbon monoxide (CO) molecules adsorbed on Cu(lll). We show how chemically 
functionalized STM tips modify the lETS intensity corresponding to adsorbate modes on the sample 
side. The underlying propensity rules are explained using symmetry considerations for both the 
vibrational modes and the molecular orbit als of the tip and sample. This suggests that single- 
molecule lETS can be optimized by selecting the appropriate tip orbital symmetry. 



I. INTRODUCTION 

Vibrational spectroscopy with the scanning tunnel- 
ing microscope (STM) is a powerful technique to study 
molecules adsorbed on surfaces.^ The method combines 
the high spatial resolution of the STM with the chemi- 
cal sensitivity of inelastic electron tunneling spectroscopy 
(JETS) as first reported for molecules incorporated in 
metal-oxide- metal tunneling junctions.^ At low tempera- 
tures, molecular vibrations are almost completely frozen 
but can be excited by tunneling electrons with sufficient 
energy. By recording the changes in tunneling conduc- 
tance as a function of applied bias voltage the onsets of 
distinct phonon emission processes can be detected. In 
this way spectra recorded on single molecules are vibra- 
tional fingerprints that reveal information about chemical 
species and local environment. Inelastic tunneling elec- 
trons from the STM can also be used to induce and con- 
trol molecular motion^ — and single-molecule chemistry^ 

Several theoretical methods have been developed to de- 
scribe lETS-STM,^^^ including approaches based on ah 
initio calculations for many-electron systems J^ — Good 
quantitative agreement between simulated and experi- 
mental data have been reported for many systems^ — 
but the spectra are difficult to rationalize due to the 
complexity of the underlying symmetries of electronic 
structure and vibrational modes. Therefore, in order to 
identify the active vibrational modes in lETS and to un- 
derstand the corresponding magnitude and sign of the 
conductance changes various theories of inelastic propen- 
sity rules have been developed^^— The many-body ex- 
tension of the Tersoff-Hamann approach developed by 
Lorente and Persson^^ permits the derivation of selection 
rules based on the symmetry of the vibrational modes 
and the molecular orbitals close to the Fermi level^ 
However, the morphology and chemical composition of 
the STM tip is often largely unknown in the experiments 
and, therefore, tip effects can be important for the inter- 



pretation of experimental data^^ — For example, chemi- 
cally functionalized tips were shown to reveal vibrational 
modes which are not observed with a bare metallic tip^^ 
Consequently, understanding and disentangling the role 
of the STM tip in lETS measurement are of relevance. 

In this paper we investigate theoretically the interplay 
between symmetries imposed by the electronic structure 
of the STM tip and the lETS spectra originating from a 
vibrating molecule adsorbed on the metal sample. Specif- 
ically, we consider the prototype system of a CO molecule 
adsorbed on Cu(lll) in combination with bare and chem- 
ically functionalized STM tips. This choice is motivated 
by the fact that this system has been the subject of re- 
search with a variety of surface analysis techniques for 
decades and that the observed lETS with metallic tips 
is now well understood ] ^^i^^i^^ We show that chemically 
functionalized STM tips can increase the resolving power 
of lETS and suggest the importance of the symmetries 
imposed by the electronic structure of the tip probe. 

The outline of the paper is as follows. In Sec. [Ill we 
describe the methodology used in the present study to 
simulate lETS. In Sec. [Till we report our calculations for 
several specific junctions. In a first set of calculations 
(Sec. nil Ap we examine aligned junctions in which the 
sample CO and the molecule on the tip side share a com- 
mon rotation axis. These results are analyzed in terms of 
the symmetry and spatial distribution of the most trans- 
mitting eigenchannels of the junction in order to estab- 

in Sec. 



lish the inelastic propensity rules. Then, in Sec. IIIIB( 
we consider junctions in which the tip has been later- 
ally displaced. This relaxes the symmetry constraints 
and, therefore, increases the number of allowed transi- 
tions contributing to the lETS signal. We find that the 
symmetry of the scattering states can still be used to 
understand the lETS signal. However, rather than the 
global symmetry of these, it is the local symmetry of the 
tip scattering states in the region where the deformation 
potential is significant that determines whether or not 
they are involved in the inelastic scattering. Finally, in 
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FIG. 1. STM setup: CO molecule on Cu(lll) with (a) a 
Cu adatom tip, or with a tip functionalized with (b) CO, (c) 
C2H4, or (d) CO 2. Throughout this paper the lower surface 
with the CO molecule is considered the STM sample side and 
the upper surface the tip side. 



Sec. II VI we summarize our study and formulate some gen- 
eral predictions for the influence of the tip symmetry in 
lETS. 



II. THEORY 

We use density functional theory (DFT) combined 
with non-equilibrium Green's function (NEGF) methods 
to perform lETS-STM simulations. The generic systems 
considered for the simulations are shown in Fig. [H We 
use a supercell description of the scattering region with a 
3x3 representation of a slab containing 6 Gu(lll) layers. 
To represent the STM sample side, we place a single CO 
molecule adsorbed on a top-site of one of the Cu(lll) 
surfaces. On the opposite surface, we place a single Cu 
adatom [Fig. [U^a)], a carbon monoxide [CO, Fig. [Hb)], 
an ethylene [C2H4, Fig. [TJc)] or a carbon dioxide [CO2, 
Fig. md)] molecule to model the STM tips. 

The electronic structure and geometries are calculated 
with the Siesta code.^^ Atomic coordinates of the CO 
molecule, tip, and surface copper atoms were relaxed 
until forces were smaller than 0.02 eV/A. The basis set 
consisted of single-zeta plus polarization orbitals for Cu 
atoms and split- valence double-zeta plus polarization or- 
bitals for C, O, and H atoms. We used the general- 
ized gradient approximation for exchange-correlation^ 
a cut-off of 200 Ry for the real-space grid integrations, 
and the F-point approximation for sampling of the three- 
dimensional Brillouin zone. 

The transport properties are simulated with the 
TranSiesta setup^ in which the scattering region de- 
scribed above is coupled to two semi-infinite electrodes 
on the sample (S) and tip (T) sides, respectively. Since 
we are only interested in the low-bias regime it suffices to 
calculate the electronic structure in equilibrium. To keep 
the present analysis of the inelastic scattering rates trans- 
parent we also restrict the transport calculations to the F- 
point of the two-dimensional Brillouin zone. Eigenchan- 
nel visualizations are calculated according to the method 
described in Ref. 34. 

To simulate the inelastic corrections to the tunneling 



current we employ the Inelastica package. ^^'^^ The vi- 
brational modes A and energies hujx are obtained by di- 
agonalizing the dynamical matrix extracted from finite 
differences for atoms in a predefined vibrational region. 
Here, we consider two cases: (i) a realistic situation where 
the sample CO molecule, the tip structure, and the Cu 
atoms in the topmost layers of Cu(lll) surfaces are con- 
sidered as vibrationally active components of the system, 
or (ii) an imaginative but illustrative situation where 
only the sample side (CO and topmost Cu layer) vibrates 
and the tip side atoms are frozen. To each vibrational 
mode we determine the corresponding electron-phonon 
(e-ph) coupling matrix M'^ via finite differences of self- 
consistent Kohn-Sham Hamiltonians^^ 

The inelastic scattering is addressed using the NEGF 
formalism within the lowest-order-expansion (LOE) 
scheme^— which relies on weak e-ph interactions and 
the electronic structure at the Fermi level only. Ignoring 
quasi-elastic correction term a^^i^^ we write the inelastic 
current as: 

/LOE = Go Vr 

+ E <^A))Tr [Gtr^ G{M^ At 

A 

+ ^(rTG'^M^ AM^ -H.c.)}] 
+ E '^^)Tr [G^Ts G 

A 

x{rTG'^M^(AT - A5)M^ + H.c.}]. (1) 

The first term in Eq. ([T]) is the elastic Landauer re- 
sult with transmission r times the conductance quan- 
tum Go = 2e^//i. The second and third terms con- 
tain the inelastic corrections to the current from each 
vibrational mode at the corresponding threshold volt- 
ages. The functions /f ^ (F, kT, (nx)) and /f (F, kT) 
give rise to symmetric and asymmetric features in the 
differential conductance dl/dV and depend on bias volt- 
age V, fundamental temperature kT, and thermal ex- 
pectation value of the mode occupation (nx). The re- 
tarded Green's function G, the spectral function A = 
i{G — G^), as well as the electrode couplings Ts,t are 
all evaluated at the Fermi energy. The quantities 
^s,T = G Ts,T G^ are partial spectral functions defined 
such that A = A^ + At- The sum in Eq. ^ runs over 
all modes A in the vibrational region. 

From the inelastic current I{V) the lETS, defined as 



lETS 



^ dl/dV 



(2) 



can be calculated. The normalization by dl/dV is con- 
venient to eliminate the exponential dependence on the 
tip-substrate distance of the tunneling current. In this 
way the normalized lETS is a measure of the inelastic 
cross sections of the various modes. Examination of the 
"asymmetric" term in Eq. ([T]) shows that it is zero for 
symmetric systems. Experimental lETS curves do of- 
ten contain asymmetric signals, although the size of these 
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signals is usually small. Such inelastic features can be 
removed by averaging the lETS signals for positive and 
negative bias. This symmetrization procedure has been 
applied to all the lETS spectra reported in Sec. Illli Fur- 
thermore, our spectra are broadened according to a lock- 
in signal Vrms = 1-8 mV and a temperature of T = 6 
We note that while the bias inversion symmetry is 
only approximate for asymmetric junctions, this symme- 
try is unrelated to the orbital and vibrational symmetries 
discussed below for the lETS propensity rules. 

Our results are useful to understand the effect of the 
symmetry on the changes in the strength of the lETS 
signals. With the aim of elucidating the origin of these 
changes, we have thoroughly investigated the role of the 
tip electronic structure on the inelastic scattering asso- 
ciated with the vibrational modes of the CO molecule 
on the sample side. Nevertheless, caution should be ex- 
ercised in the quantitative comparison of our simulated 
lETS spectra with experiments, since various approxi- 
mations are involved in our calculations, such as F-point 
sampling and LOE approach. Although for most cases 
our approximations are reasonable, in some situations 
they might affect position and amplitude of the lETS 
peaks. 



III. RESULTS 

In our calculations we have considered the different 
junctions shown in Fig. [H The sample is always a flat 
Cu(lll) surface with an adsorbed CO molecule. The 
molecule is assumed to adsorb vertically on a top-site 
of Cu(lll), as experimentally observed. Although such 
assumption overcomes some of the limitations of many 
DFT functionals which incorrectly give the fcc-hollow as 
the most stable configuration,—^ it introduces inaccura- 
cies in the prediction of the low-energy frustrated trans- 
lation (FT) mode.^^ Therefore, in the following we will 
focus on the frustrated rotation (FR) and other higher- 
lying vibrational modes, well-defined in our DFT-NEGF 
calculations/^ 

Three conductance channels dominate the elastic 
transport through the CO molecule with different rota- 
tional symmetry around the molecular axis. As shown 
in Fig. [21 two channels have tt character and one has 
a character^^ Since electron propagation can occur 
through each channel in either direction (from sample 
to tip or vice versa) there are thus two different current- 
carrying scattering states or wave functions (denoted Si 
and ti) associated with the z'th channel. In the STM 
setup these scattering states are exponentially damped 
across the tunneling barrier and retain only a tiny ampli- 
tude on the transmitted side (related to the transmission 
coefficient). Among the three dominant sample states 
{si, ^2 , S3 } (i.e., CO side of the junction) we see in Fig. [2] 
that the weight of the state is strongly reduced in the 
molecular apex (O atom) compared with the more ex- 
tended states. As expected, the symmetry and spatial 



distribution of the sample states reflect to some extent 
those of the frontier orbit als of CO: the two- fold degener- 
ate TT-type LUMO, and the a-type HOMO which is more 
localized on the C atom. 

The symmetry of the channels is an important ingre- 
dient to determine the inelastic cross-section associated 
with a given vibrational modej^^ Given the symmetry 
of the FR mode (tt character with respect to the CO 
molecule symmetry axis), cf. Fig. [2fe), we can formulate 
the selection rules in a quite concise form: sample states 
with TT character couple only to tip states with a char- 
acter around the axis of the sample CO molecule (i.e., 
^ tj) and, conversely, sample states with a char- 
acter couple only to tip states with tt character (i.e., 
5^ ^ t^)M Note that here we always define the tt or 
a symmetry character of the transmission channels with 
respect to the axis normal to the surface, i.e., with re- 
spect to the molecular axis, and not in the conventional 
way for molecular orbit als. Therefore, its meaning might 
be different, in particular when the two symmetry axes 
do not coincide. 

It is important to note that symmetry considerations 
only allow to determine which scattering events are pos- 
sible. The relative contribution of each allowed transi- 
tion to the total inelastic cross-section depends on other 
factors. In particular, inelastic scattering rates strongly 
depend on the extension of the wave functions of the 
scattering states (both from sample and tip) along the 
molecule and towards the vacuum. Since we focus on 
molecular vibrations, the cross section of each scattering 
event is determined by the overlap between these scat- 
tering states within the gap region, where the molecule 
is located and the deformation potential associated with 
the vibration is stronger. From this simple reasoning, 
and given their different spatial distribution, we can ex- 
pect that transitions in which of the CO molecule is 
involved will be less efficient that those involving chan- 
nels of the same molecule. As we will see, this is indeed 
confirmed by our calculations. 

In the present theoretical work we analyze the ef- 
fect of changing the symmetry and extension of the tip 
states via functionalization of a metallic tip with different 
molecules, cf. U states in Fig. O Our model tip is com- 
posed by a flat Cu(lll) surface with or without a Cu 
adatom, or with one of these three different molecules: 
CO, CO2, and C2H4. In a first set of calculations, the tip 
apex is always aligned with the sample CO molecule, so 
both sample and tip share a common symmetry axis. The 
use of these aligned sample-tip arrangements allows us to 
investigate the effect of the intrinsic symmetry of the tip 
states on the inelastic transport properties using selection 
rules. In a second set of calculations, the tips are laterally 
displaced. The main effect of this displacement is to relax 
the symmetry constraints of the system, increasing the 
number of allowed transitions, and to change the distance 
between the sample and tip molecules, thus affecting the 
overlap between scattering states coming from both sides. 
Although the lateral displacement strongly affects the in- 
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FIG. 2. Isosurface plots of the three most transmitting eigenchannels for each of the four tips considered: 
(a) Cu adatom tip, (b) CO tip, (c) C2H4 tip, and (d) CO2 tip. The states labeled Si (ti) correspond to 
electron waves incoming from sample/left (tip/right) side. Due to the tunnel gap between tip and sample 
the electron scattering states decay rapidly and the transmitted part of the wave on the other electrode is 
not visible. The blue and white colors represent the sign of the real part of scattering states (our choice 
of phase makes the imaginary part negligible for visualization purposes), (e) Visualization of the FR 
vibrational mode of CO. 



elastic transport in some cases, it is important to note 
that such changes can still be completely rationalized us- 
ing the same ingredients that allow to understand the 
more symmetric situations: (i) the local character (a or 
tt) of the most conducting tip states around the sam- 
ple CO axis, and (ii) the overlap between sample and 
tip states with the deformation potential. All the re- 
sults reported in the following correspond to the tunnel- 
ing regime, with distances between the most protruding 
tip and sample atoms > 3 A. 



A. Aligned tips 



Inelastic spectra for CO/Cu(lll) using a Cu adatom 
tip have recently been experimentally and theoretically 
characterized We will use this system as a standard 
reference in what follows. The generic lETS spectrum is 
dominated by the FR of the CO molecule at ~ 35 meV, 
cf. Fig. [3l At lower energies, not considered here, the 
FT of CO gives a weaker signal at 5 meV. Similarly, in 
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FIG. 3. Symmetrized lETS spectra of a CO molecule ad- 
sorbed on a Cu(lll) surface with different models for the 
tip: (i) a CO functionalized tip that is either vibrating (solid 
line) or frozen (dashed line), (ii) a vibrating flat Cu(lll) tip 
(dashed-dotted line) or (iii) a vibrating sharp Cu adatom tip 
(dotted line). As a reference, the calculated relative change 
in tunneling conductance due to the frustrated rotation (FR) 
modes is 24% with the Cu adatom tip but only 3% with the 
flat Cu tip. The inset shows a zoom on the signal related to 
the CO center-of-mass motion around 57 mV. 



our simulations with functionalized tips presented below 
we also find the spectra dominated by the FR of the CO 
molecule at ~ 35 meV (see Fig. 2]) for all the different 
tip terminations, albeit with varying intensity. However, 
there are also important differences that deserve a careful 
analysis. 

First, we examine the case of the CO functionalized 
tip, where the top- most Cu(lll) layers and both tip and 
sample CO molecules are allowed to vibrate. The corre- 
sponding theoretical lETS spectrum is shown in Fig. [3] 
together with the result for the Cu adatom tip. In ad- 
dition to the FR signal, we find a weaker feature that 
originates in the center-of-mass (CM) motion (Cu-CO 
stretch) at ~ 57meV (see inset to Fig. [3]). Clearly, the 
use of a CO functionalized tip enhances the CM signal, 
which was hardly observed with a Cu adatom tip This 
effect follows from the selection rules described in Ref. 
and the observation of such an enhancement in the CM 
signal might be employed for the characterization of the 
STM tip used in experiments. 

In order to check if the presence of the two CO 
molecules at the relatively short tunneling distance used 
here introduces some coupling of the vibrations or non- 
linear effect in the electron-phonon coupling, we have 
performed a calculation in which one side of the junc- 
tion is kept frozen. This asymmetric vibrational region 
yields asymmetric features in the lETS spectrum which 
we here average out as described above. As seen in Fig.[3l 
the FR signal obtained when both tip and sample CO 
molecules vibrate (black line) is larger than the signal 
obtained with a frozen CO tip (blue dashed curve) by a 



factor of ^ 2, showing that the coupHng between the vi- 
brational modes of the two molecules, and its effect in the 
e-ph interaction, is negligible at the considered distances. 
Therefore, the deformation potential associated with the 
CO FR mode is not substantially affected by the pres- 
ence of neighboring vibrating molecules, and the changes 
in the calculated lETS spectra using different function- 
alized tips are mainly related to the different symmetry 
and spatial distribution of the electronic states in each 
tip. 

Using a CO functionalized tip and a Ag(llO) substrate, 
Hahn and Hc^^ measured the relative intensity of the FR 
signal for a CO- vacuum- CO tunnel junction and for a 
CO-vacuum-Ag(llO) tunnel junction. The inelastic sig- 
nal was observed to be considerably larger in the former 
case. In order to check if we can reproduce this behavior, 
we have also calculated the lETS spectrum using a fiat 
Cu tip represented by a Cu(lll) surface. Figure [3] shows, 
in agreement with the experiment, that the correspond- 
ing FR signal (red dash-dotted curve) is smaller than the 
signal obtained with a CO tip (black line). We note that 
the calculated signal at low energies (around ±20 mV) 
originates from the vibrational modes of the Cu surface 
and its analysis is out of the scope of this work4^ 

Next we consider the two other tip terminations shown 
in Fig. [H As opposed to the CO molecule which is ori- 
ented perpendicular to the Cu(lll) surface, the C2H4 
and CO2 molecules form almost planar configurations 
on Cu(lll)4^— The corresponding lETS spectra are 
represented in Fig. IH where the positive voltage region 
(V > 0) corresponds to a realistic vibrational region 
which includes the sample CO, the tip and the top- most 
Cu layers, while the negative voltage region (V < 0) 
shows the results obtained for the imaginative but illus- 
trative case when only the sample side (CO molecule and 
topmost Cu layer) vibrates. A comparison of the signals 
shown for these two types of vibrational regions allows us 
to conclude that the lETS signal around ±35 mV origi- 
nates from the FR modes of CO for all four tip models 
(the signal amplitude is practically unchanged by freez- 
ing the tip atoms). Furthermore, with the CO2 tip frozen 
the CO CM signal shown around — 60mV is suppressed 
(see inset to Fig. |4]). The larger signal observed with the 
CO2 tip vibrating (shown for V > 0) thus originates in 
tip vibrations. The FR signal is stronger for the CO2 tip 
than for the Cu adatom tip, which in turn is much larger 
than the signal with CO or C2H4 tips. A small FR sig- 
nal with the C2H4 tip is consistent with the experimental 
observations reported in Ref. [28|. 

Now we discuss our findings using the above mentioned 
symmetry considerations for both the vibrational mode, 
more precisely the corresponding deformation potential, 
and wave functions of the relevant (tip and sample) scat- 
tering states, as well as their spatial extension. The in- 
elastic scattering rate for a given phonon mode A can be 
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FIG. 4. Symmetrized lETS spectra of a CO molecule ad- 
sorbed on Cu(lll) surface using a Cu adatom tip or a tip func- 
tionalized with CO, C2H4, or CO 2. Two vibrational regions 
are considered: (i) top-most Cu layers, sample CO molecule 
and tip molecule vibrating (shown for V > 0), and (ii) only 
sample side vibrating (shown for < 0). As a reference, 
the calculated relative change in tunneling conductance due 
to the frustrated rotation (FR) modes is 25% with the CO2 
tip but only 2% with the C2H4 tip. The inset shows a zoom 
on the signal related to the CO center-of-mass motion around 
-57 mV. 
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where represents the e-ph coupling and the sum- 
mation runs over all sample (si) and tip (t^) scattering 
states. As reported in Ref. in this eigenchannel basis 
it suffices to include only the most transmitting eigen- 
channels for the evaluation of Eq. (j3j) . This is due to the 
fact that only these scattering states have significant spa- 
tial weight in the molecular part of the junction where 
the e-ph interaction (i.e., deformation potential) associ- 
ated with the molecular vibrations is localized. In Fig. [2] 
we visualize the three most transmitting eigenchannels as 
represented by scattering states originating from the CO 
sample side {si) or from the tip side (ti). By construction 
the scattering states belonging to a given eigenchannel 
have the same symmetry. As one can see in Fig. [2l the 
three dominant channels are either a- or 7r-type. 

For each vibrational mode it follows from Eq. ([3]) that 
only scattering states with certain symmetry can couple, 
giving rise to the so-called selection or propensity rules. 
The two-fold degenerate FR mode corresponds to the ro- 
tation of the CO molecule parallel to the xy plane. The 
associated deformation potential has a tt character with 
respect to the CO axis. Correspondingly, the FR mode 
can only produce scattering of electrons between tt- and 
cr- type channels. This simple picture is confirmed by our 
numerical calculations for the aligned junctions. Table |I] 
gives the calculated contributions of the transitions be- 
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TABLE I. Decomposition (in %) of the total inelastic scat- 
tering rate for the frustrated rotations (FRi and FR2) of CO 
on Cu(lll) among the three dominant scattering states only 
(i.e., a 3 X 3 matrix). The sum over all 9 matrix elements, 
stated next to the tip label, does not reach 100 % because 
of the total inelastic scattering rate also involves scattering 
with the other, less-conducting channels not analyzed in this 
work. The elastic transmission coefficient for each channel is 
also given (ji — — rl). 



tween the three most transmitting eigenchannels to the 
total inelastic scattering rate of the two FR modes of CO, 
as obtained with different tips. For the Cu adatom tip the 
main contribution to the inelastic scattering comes from 
the transitions {sj,^^} ^ ^3 and amounts to ~ 95% of 
the total scattering rate. Similarly, for the CO2 tip the 
main contribution is {sj , S3 } ^ t J (^ 98%). In these two 
cases the dominant tip state has <j character, while the 
sample states have tt character. On the contrary, in the 
case of the C2H4 tip the dominant transition is S2 ^ t\ 
(~ 73%). For the CO tip, since sample and tip subsys- 
tems are identical, the dominant contribution (~ 94%) 
is given by the sum of transitions between {sj , } ^ ^3 
and between S3 ^ {^1,^2}- 

In the case of the Cu adatom and CO tips [see 
Fig. [Ha,b)] there exists an axial symmetry in the junc- 
tion due to their vertical adsorption and the hexagonal 
symmetry of the substrate. This implies that the two tt 
eigenchannels and the two FR modes are identical. For 
the Cu adatom tip, the tip state t%, extends more to- 
wards vacuum than and tj? cf. Fig.[2fa). In principle, 
we can relate the a channel to the 45 orbital of the Cu 
adatom, for which we can expect to find a considerable 
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spectral weight at energies nearby the Fermi energy. Ac- 
cording to the selection rules, this tip state can couple 
with the sample states {sj, }. Since and also ex- 
tend the most along the molecule and into vacuum, the 
transitions {sj,^!} ^ (~ 95%) completely dominate 
the inelastic scattering. As expected by symmetry, each 
of these two transitions contribute almost equally to the 
total inelastic rate. For the same reason, the transitions 
53 ^ {tj , ^2 } 2%) also contribute equally, although in 
a much smaller percentage due to the smaller overlap of 
S3 with the vibrational region and with {^1,^2} of 
adatom tip. With the same axial symmetry, the situation 
is similar in the case of the CO tip, where the transitions 
S3 ^ tl, si ^ ^ tl, and si ^ q yield a - 23% 

contribution each, i.e., in this case tip and sample states 
with both symmetries become relevant. However, since 
with two CO in the junction the a-type states are always 
rather localized near the metallic surface, the inelastic 
signal is significantly reduced with respect to that of the 
Cu adatom tip. This is so, even if we take into account 
that, with two vibrating CO molecules in the junction, 
the inelastic signal associated with the FR modes of CO 
doubles with respect to that of a single vibrating CO 
molecule. 

Although in this paper we mainly focus on the FR 
mode of CO, the reason for the intensity increase of the 
CM mode when a CO tip is used also deserves some dis- 
cussion. Contrary to the case of the FR mode, the de- 
formation potential associated with the CM mode has a 
character around the molecular axis and induces transi- 
tions between tip and sample channels of the same sym- 
metry, i.e., sf ^ or s^ ^ tj. Using a CO tip, the 
inelastic scattering induced by the CM mode is domi- 
nated by the more extended tt channels of CO. Using a 
Cu adatom tip, the most important transition induced 
by the CM mode is S3 o t^. However, these a chan- 
nels have a much smaller overlap as compared with the tt 
channels that dominate the case of the CO tip. For this 
reason, the CO tip considerably enhances the CM signal. 

As for CO2 and C2H4 tips, they exhibit in-plane ad- 
sorption on Cu(lll) and hence, define a preferential di- 
rection given by the molecular axis within the xy plane. 
Consequently, the two tt channels are not identical any- 
more and the scattering cross-section of the two FR 
modes also differs. In the case of the CO2 tip the most 
important processes are the transitions {^2,53} ^ 
(~ 98%), i.e., between sample tt states and tip a states, 
just as in the case of the Cu adatom tip. However, the 
contribution from 53 ^ t J is larger than the contribution 
from S2 ^ ti. This symmetry breaking is also observed 
for the Si ^ component, which is almost negligible 
in comparison with the s^ ^ component. In the case 
of C2H4, we find that the transition s^ ^ tj amounts 
to ^ 73% of the total inelastic rate while s^ ^ q only 
gives ~ 5%. This difference is readily understood from 
Fig. [2fc) where one notes that q extends more towards 
the CO sample than the q state, and hence has the larger 
overlap with the vibrational region. 



All these results provide numerical evidence of the se- 
lection rules governing the inelastic scattering processes. 
Furthermore, as the e-ph coupling is local in space, even 
if the symmetry of the orbitals allows for a transition, the 
sample and tip states must also extend quite appreciably 
in the transport direction in order to yield a significant 
inelastic signal. Therefore, we have seen that it is the 
combination of the symmetry and spatial extension of 
the tip and sample states that ultimately determines the 
lETS intensity. This is indeed what is reported in Fig. 21 
where the FR signal is smallest for the C2H4 tip, due to 
a strong attenuation of the state along the molecular 
axis as compared to, e.g., for a Cu tip. On the con- 
trary, for the CO2 tip the state extends deeper into the 
vibrational region and thus allows for a stronger inelastic 
FR signal. 



B. Laterally displaced tips 

In the following we consider junctions with laterally 
displaced tips, shown in the lower panels in Fig. [5l If the 
tip is slightly displaced from the aligned arrangements 
considered in Sec. IIII Al [see Fig. [1], the symmetry of the 
tip-sample system is broken. For the Cu adatom or CO 
tips this means that we no longer have a common rota- 
tion axis in these junction. Similarly, for the CO2 and 
C2H4 tips (which do not carry rotational symmetry) it 
means that we reduce the number of symmetry opera- 
tions. As illustrated in Fig. [5l we are considering the fol- 
lowing laterally displaced tip-sample configurations: (i) 
hollow and next-top configurations for atomic Cu, CO 
and CO2 tips, and (ii) short-bridge and long-bridge hol- 
low configurations for the C2H4 tip. We note that the 
Cu(lll) surface and the tip cannot easily be laterally 
displaced due to the use of periodic boundary conditions 
in Siesta. This would require the use of a substantially 
larger unit cell (e.g., the use of two slabs instead of one, or 
the inclusion of a dislocation in a sufficiently thick slab, 
such that the dislocation does not affect the surface elec- 
tronic structure). The use of such a large supercell would 
significantly increase the computational cost of the cal- 
culations. Instead, we have decided here to use a single 
metallic slab. This constraint keeps the computational 
load at a moderate level and facilitates the comparison 
with calculations of the aligned configuration, although 
forces us to study only discrete translations of the func- 
tional part of the tip from one adsorption site to the 
nearest one. 

The results for the corresponding inelastic spectra are 
also shown in Fig. [5l One sees that the intensity of the 
inelastic signal strongly depends on the position of the 
tip. In some cases, we also see slight variations of the 
energy of the CO FR mode. These small variations are 
consistent with those recently observed in Ref. [l8| when 
the size of the tunneling gap is changed. However, in the 
following we concentrate only in the reasons behind the 
change of the intensity of the inelastic signal. 
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FIG. 5. Symmetrized lETS spectra for a CO molecule adsorbed on Cu(lll) surface using (a) a Cu adatom tip, (b) a CO, (c) a 
C2H4, or (d) a CO2 tip placed on different positions of Cu(lll), as shown at the bottom part of each figure. Two vibrational 
regions are considered: (i) top- most Cu layers, CO molecule and tip molecule (shown for V > 0), and (ii) only sample CO 
molecule (shown for V < 0). As a reference, the calculated relative change in tunneling conductance due to the frustrated 
rotation (FR) modes is 24% with the Cu adatom tip (top site). The (green) open circle indicates the position of the sample 
CO molecule, which always remains on a Cu top site in the opposite surface. 



In the case of the Cu adatom and the CO2 tips [see 
Fig. [5ja,d)], the corresponding lETS spectra exhibit al- 
most no change when displacing the tip from top to hol- 
low position. As already shown in Sec. IIII A\ the most 
important tip states for FR scattering with these two 
types of tips are the states with a symmetry. There- 
fore, as the tip moves away from the top position, the 
main effect is a modification of the overlap between a t^- 
type state and the deformation potential associated with 
the FR modes. Such an overlap does not change appre- 
ciably as the Cu adatom or CO2 tips move from top to 
hollow positions and, as a result, the corresponding FR 
signal is hardly modified. However, when the tips are 
further displaced to the next-top position the overlap de- 
creases much more and, consequently, the FR signal is 
reduced. In principle, due to the reduction of symme- 



try the <j and tt classification is only approximate and 
all transitions among the scattering states {si ^ tj) are 
allowed. But these contributions are small in practice. 
Our calculations suggest that when a tip state gives 
the largest contribution to the inelastic signal of the FR 
mode of CO, this signal will decrease as the tip is laterally 
displaced following the reduction of the overlap between 
the electronic states of the tip and sample. 

The lETS spectra with CO and C2H4 tips present 
a different behavior as the tip-sample relative positions 
change [see Fig. [5fb,c)]. The main difference, as com- 
pared to the previously discussed cases of Cu adatom 
and CO2 tips, is that, the FR signal increases when the 
tip is laterally displaced away from the aligned position. 
The main factor behind this behavior is the increase in 
the number of allowed transitions when the symmetry is 
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reduced. In the case of CO, as the tip is laterahy dis- 
placed from its original top position to hollow site and 
further to next top, the FR signal first increases signifi- 
cantly and then decreases practically to its initial value. 
The first increase is easy to interpret: when the molecule 
is displaced inelastic scattering between tip and sample 
TT-like states becomes allowed (tt o tt). Since the tt states 
extend the most into the tunneling gap this transition is 
more efficient than a ^ tt, which were the only allowed 
transitions for the aligned configuration. Besides, there 
is a signal at ~ ±20 mV caused by vibrations of the Cu 
atoms in the Cu(lll) surface layers of the sample and the 
tip, similar to the one observed with a flat Cu(lll) tip 
[see Fig. [3]. By further displacing the tip to the next top 
position we see a decrease in the FR signal, consistent 
with a decrease of the overlap between tip states with 
sample states as well as the vibrational region. 

A similar behavior is observed in Fig. [5fc) in the case 
of C2H4 functionalization. When the tip molecule is lo- 
cated on top site or long-bridge hollow site the FR signal 
is almost identical and small, but dramatically enhanced 
if we move the tip to a short-bridge hollow position. In 
order to understand the origin of this strong enhance- 
ment of the FR peak we analyze the spatial properties 
of the dominant tip state {ti) for each of the three C2H4 
tip positions. Figure [6] shows a top view of the three 
different tip states, where the origin (green empty circle) 
corresponds to the position of the sample CO molecule. 
This allows us to survey the "local" symmetry character 
of the tip states with respect to the center of the CO vi- 
brational mode. As mentioned in the beginning of Sec. lIIII 
this is the critical feature responsible for the amplitude 
of the lETS signal of the FR mode. In Figs. [6l^a,c) we 
observe that for top and long-bridge hollow tip arrange- 
ments, the dominant tip state has a local 7r-type symme- 
try. For this reason the FR mode cannot couple it effi- 
ciently to the TT-type CO sample states, hence the small 
FR inelastic signal in Fig. [5fc). On the contrary, when 
the C2H4 molecule is placed in the short-bridge hollow 
position [see Fig. [6fb)] the dominant eigenchannel has a 
strong local a-type character and therefore yields a sig- 
nificantly stronger FR signal. The lETS peak for the FR 
modes of the CO is even larger than the one observed 
with a Cu adatom or a CO2 tip on top positions, due to 
the fact that the involved tip states extend even further 
into the vibrational region. 



IV. CONCLUSIONS 

In this work the role of the STM tip on inelastic elec- 
tron tunneling spectroscopy has been thoroughly inves- 
tigated. Special emphasis has been placed on the effects 
of changes in symmetry of the orbitals involved in the 
inelastic scattering processes. 

We have simulated lETS spectra for a CO molecule 
placed on a Cu(lll) surface using a Cu adatom tip as 
well as tips functionalized with different molecules (CO, 




FIG. 6. Top view of the tip scattering state which contributes 
the most to the inelastic scattering rate. The green empty 
circle represents the xy position of the sample CO molecule. 
The panels correspond to the C2H4 molecule placed on (a) 
top site, (b) short-bridge hollow, and (c) long-bridge hollow 
positions [as shown in the lower panel of Fig. [SJ^c)]. For clar- 
ity, the Cu atoms on the first layer underneath the tip have 
been represented in dark brown. The blue and white colors 
represent the signs of the real part of scattering states (due 
to the choice of phase the imaginary part is negligible for 
visualization purposes) . 



CO2, and C2II4). For all the different tips considered 
here, the lETS spectra are dominated by the excitation 
of the frustrated rotation vibrational mode of the CO 
molecule at '^35 meV. However, we find that the intensity 
of the FR peak strongly depends on the chemical nature 
as well as on the position of the STM tip. These effects 
were rationalized by means of the propensity rules, i.e., 
taking into account the local symmetry character of both 
the deformation potential associated with the excited vi- 
brational mode and of the most transmitting scattering 
states associated with the tip and sample, in the region 
where the deformation potential is appreciable. 

We performed a detailed analysis of these effects by 
visualizing the eigenchannels of the different junctions. 
In most cases the symmetry and the spatial distribution 
of these dominant scattering states can be rationalized 
in simple terms and correlated with the frontier molec- 
ular orbitals. Specifically, our analysis suggests the fol- 
lowing consequences of lETS propensity rules related to 
the symmetry respect to the surface normal: Adsorbate 
molecules for which 7r-type frontier orbitals control the 
electron transport reveal lETS signals from transversal 
(longitudinal) vibrational modes if observed with tips of 
cr-type (TT-type) states. Conversely, adsorbate molecules 
with dominating cr-type frontier orbitals reveal lETS sig- 
nals from longitudinal modes if observed with cr-tips and 
transverse modes if observed with 7r-tips. In short, ro- 
tational symmetry around the transport direction favors 
inelastic transitions where the product of sample state, 
deformation potential, and tip state is symmetric. Note 
that the tt or a symmetry character of the transmission 
channels is defined with respect to the axis normal to 
the surface and, therefore, its meaning is not always the 
same as the conventional one of the molecular orbitals. 

Our calculations further show that, for a given vibra- 
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tional mode of the adsorbed molecule, the inelastic scat- 
tering rate is mainly determined by the local character 
in the vibrational region of the tip channels and that, 
as a result, the relative tip-sample configuration is cru- 
cial. Therefore, in the general cases where sample and 
tip do not share a symmetry axis, we have shown that 
the propensity rules stated above apply approximately 
in terms of the "local" symmetry in the region of space 
where the molecular vibrations occur. 

For the particular case of the frustrated rotation of CO, 
the locality effect is specially relevant if the tip states 
have dominant ly 7r-type character, because then the in- 
elastic signal can be largely increased by simply shifting 
the tip laterally. We have also found that, in some cases, 
functionalized tips can enhance inelastic signals as com- 
pared to those obtained with a metallic adatom tip. For 
instance, the signal from the center-of-mass motion of 
one CO molecule was found to be larger when observed 
with a CO tip than with a metallic tip, cf. inset to Fig.[3l 

In conclusion, we have shown that the STM tip is of 



vital importance in inelastic tunneling spectroscopy and 
that, as a result, single-molecule vibrational spectroscopy 
could be modulated and manipulated by means of func- 
tionalized STM tips. Moreover, since the effect of the tip 
might play an important role in the JETS measurements, 
it would be desirable to take it into account in the inter- 
pretation of experimental data. We believe our results 
can help in that direction. 
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